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Processing Seismic Data 



Field of the Invention 

5 The present invention relates to a method of processing seismic data, in particular to 
identify seismic events arising from microseismic activity, for example microseismic 
activity in a hydrocarbon reservoir, a nuclear waste repository, or a mine, or in 
seismology. 



10 Background to the Invention 

Figure 1 is a schematic view of a hydrocarbon reservoir. In essence, a hydrocarbon 
reservoir consists of layers of the earth's interior in which hydrocarbons have 
accumulated. Only one hydrocarbon-bearing layer 2 is shown in Figure 1, for 
15 convenience of description. In Figure 1 the structure of the earth overlying the 
hydrocarbon-bearing layer 2 is represented as a single layer 1, but in practice it will 
consist of many layers having different geological compositions. Similarly, the 
structure of the earth below the hydrocarbon-bearing layer 2 will in practice consist of 
many layers but is represented in Figure 1 as a single layer 3. 

20 

Microseismic activity, or "microseismicity", in hydrocarbon reservoirs can be induced 
as a result of extraction of hydrocarbons from the reservoir, which has the effect of 
altering the weight distribution of the hydrocarbon-bearing layer 2 and of the overlying 
layer 1. Microseismicity may also be induced by hydraulic fracturing operations. 

25 When microseismic activity occurs, seismic waves are generated at a point within the 
earth's interior. The seismic waves propagate through the earth, and may be detected by 
seismic receivers located within the earth. Microseismic activity induced by 
hydrocarbon extraction is weak and can not generally be detected by seismic receivers 
at the earth's surface, although earthquake-induced microseismic activity can generally 

30 be detected by receivers at the earth's surface. 



In principle, microseismic activity can occur anywhere in the reservoir layer 2 of Figure 
h in the overlying layer 1 or in the surrounding rock. In practice, however, the earth's 
structure contains geological fractures or other faults, and microseismic activity 
preferentially occurs along or in the vicinity of such faults or fractures. In Figure 1, a 
fault 4 is schematically depicted in the hydrocarbon bearing layer 2. It is assumed that 
the earth's interior to the right of the fault 4 is tending to move downwards whereas the 
earth's interior to the left of the fault 4 is tending to move upwards as indicated by the 
arrows A,B in Figure 1. If a portion of the fault at one location 5 slips, microseismic 
activity occurs at that location 5. The slippage will relax tension at the location 5 on the 
fault where slippage occurs, but will induce tension at neighbouring locations along the 
fault and this increased tension may give rise to slippage, and consequent microseismic 
activity, at a nearby location 6 on the same fault 4. The time delay between 
microseismic activity at one location 5 and microseismic activity at the nearby location 
6 may be of the order of seconds or minutes, or it may be of the order of days or even 
weeks. 

The seismic energy produced by the microseismic activity occurring at locations 5,6 
will be detected by seismic receivers 7,8,9 disposed in a borehole. 10. The seismic data 
acquired at the receivers 7,8,9 therefore contain events arising from microseismic 
activity - or "microseismic events" - in addition to microseismic events arising from 
other faults in the layers 1,2. 

Where microseismic activity occurs along a geological fault or fracture, microseismic 
activity is, as explained above, often found to occur at two nearby locations, such as the 
locations 5,6 in Figure 1. Since the two locations are near to one another, the focusing 
effects of the overlying layers on seismic waves emitted at one location 5 will be similar 
to the focusing effects of the overlying layers on seismic waves emitted at the nearby 
location 6. Furthermore, since the microseismic activity at each location 5,6 arises from 
slippage of the fault 4, microseismic activity at one location 5 has the same source 
mechanism as microseismic activity at the nearby location 6. Thus, microseismic 
activity occurring at location 5 and microseismic activity occurring at the nearby 
location 6 would produce similar recordings at a seismic receiver, since the two 
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locations are close to one another and have similar focusing mechanisms to one another, 
and since the roicroseismic activity has the same source mechanism at each location. 
Data acquired at a receiver will therefore contain an event corresponding to 
microseismic activity at location 5 and an event corresponding to microseismic activity 
5 at location 6. The two events will have a similar form, and will be separated by a time 
delay of, typically, from seconds to days or even weeks. Such pairs of events in 
acquired seismic data are known as multiple acoustic emissions or "doublets". The 
effect of the relative positions of the receiver and the event locations 5,6 on the time 
delay between the two events acquired at a receiver is small, typically of the order of 
10 10" 1 to 10" 2 seconds, so the time delay is determined primarily by the time delay 
between microseismic activity at one location and microseismic activity at the other 
location. (The effect of the relative locations of the receiver and the events on the time 
delay between the two events arises primarily from a change in the azimuth of the 
receiver with respect to the doublet orientation.) 

15 

It is desirable to identify doublets in seismic data acquired at seismic receivers in a 
hydrocarbon reservoir. As is known, once a doublet has been identified it is possible toi 
determine the relative location of the microseismic activity giving rise to each event (the; 
location of the microseismic activity giving rise to an event will be referred to as the 
20 "location of an event", for convenience). The relative locations of the two events of a 
doublet - that is, the location of one event of the doublet relative to the location of the 
other event - may be determined more precisely then their absolute locations. The 
events and their accurate relative locations can be used for fault delineation, 
permeability estimation, tracking of flood fronts, stress transfer, etc. 

25 

Acknowledgement of the Prior Art 

Initially, doublets in seismic data were identified manually. This is possible in 
earthquake studies. However, during hydrocarbon extraction a continuous seismic 
. 30 record is kept, and this record may extend for months, even years - so manual 
identification of doublets is not possible. 
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Techniques for determining the relative locations of doublet acoustic emissions have 
been reported by, for example, G. Poupinet et all in "Monitoring Velocity Variations in 
the Crust Using Earthquake Doublets: An Application to the Calaveras Fault, 
California", Journal of Geophysical Research, 89, 5719-5731 (1984), by H. Moriya et 

5 al. in "Precise source location of AE doublets by spectral matrix analysis of the triaxial 
hodogram": Geophysics, 36-45 (1994), and by F.M. Waldhauser and W.L Ellsworth in 
"A Double-Difference Earthquake Location Algorithm: Method and Apphcation to the 
Northern Hayward Fault, California", Bulletin of the Seismological Society of America, 
90, 1353-1368 (2000). However, there is currently no available method that enables 

10 doublets in a large seismic dataset to be identified quickly and easily with minimal 
operator intervention. 

Statement of the Invention 

15 A first aspect of the present invention provides a method of processing seismic data, the 
seismic data comprising at least first seismic data traces acquired at a first seismic 
receiver and second seismic data traces acquired at a second receiver spatially separated 
from the first receiver, the method comprising: determining an overall measure of 
similarity for a pair of events in the seismic traces, the overall measure of similarity 

20 being indicative of similarity between the events acquired at the first seismic receiver 
and of similarity between the events acquired at the second seismic receiver. 

When microseismic activity at two nearby locations gives rise to. a doublet in acquired 
seismic data, the doublet will appear in the seismic data as two similar events with a 

25 time delay between them. The present invention makes use of the fact that, where 
events due to microseismic activity occur in seismic data acquired at two or more 
receivers, the microseismic activity will give rise to corresponding doublets in the 
seismic data acquired at each receiver. Thus, if a doublet exists there should be a high 
degree of similarity between the two events in the data acquired at the first receiver 

30 (indicative of the occurrence of two events that are similar to one another in the or each 
data trace acquired at the first receiver) and a high degree of similarity within the data 
acquired at the second receiver (indicative of the occurrence of two events that are 
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similar to one another in the or each data trace acquired at the second receiver). The 
overall measure of similarity may be simply compared with a pre-set threshold, to 
provide a simple test for the presence of a doublet. Decomposition of the acquired data 
into P- and S- constituents is not required. 

5 

The method of the present invention may be automated, and performed with minimal 
operator intervention. This allows doublets in acquired seismic data to be identified 
quickly and easily. As explained above, acquired seismic data will typically contain 
events arising from many faults, and the invention makes it possible to identify events 
10 arising from the same fault. 

A multiplet in seismic data is a group of doublets which have events in common - the 
simplest multiplet is a pair of doublets in which one event is common to both doublets. 
Once doublets have been identified by the method of the present invention, it is then 
15* possible to identify multiplets by grouping doublets that contain a common event. 

*■/ 

The present invention makes it possible to identify doublets in real-time or in near real- 
time. 

20 In a preferred embodiment the measure of similarity is a correlation coefficient 
indicative of correlation between the events acquired at the first seismic receiver and of 
correlation between the events acquired at the second seismic receiver. 

In a preferred embodiment the method comprises: 
25 (a) determining a first correlation coefficient for the pair of events from the first 
data traces; 

(b) determining a second correlation coefficient for the pair of events from the 
second data traces; and 

(c) determining the overall correlation coefficient for the pair of events from the 
30 first correlation coefficient and the second correlation coefficient. 



When two events in the first seismic data traces (i.e., in the data acquired by the first 
receiver) are similar to one another, the first correlation coefficient will generally have 
a maximum value at a certain time lag. If the two events are a doublet, then data 
acquired at a nearby receiver should also contain a pair of events that have the same, or 
similar, time lag - since the microseismic activity should give rise to corresponding 
doublets in the data acquired at each receiver. The principle of the present invention is 
to identify events occurring with the same time lag (within a time window 
compensating for the difference between the time lag at one receiver and the time lag at 
the other receiver arising from the separation between receivers) in data acquired at each 
of two or more receivers. A pair of events that are present in the seismic data acquired 
at the first receiver but that have no corresponding events in the seismic data acquired at 
the second receiver (and so are unlikely to represent a doublet) will be given low weight 
in step (c). Thus, making use of the correlation coefficient of seismic data acquired at 
one receiver and the correlation coefficient of seismic data acquired at a second receiver 
provides more accurate identification of a doublet than if data from just one receiver is 
used. 

Furthermore, the two events of a doublet normally have very similar shapes in acquired 
seismic data traces. As a result, a doublet generally has a high correlation coefficient. 
Comparing the overall correlation coefficient with a pre-determined threshold, and 
rejecting pairs of events having a cross-correlation coefficient below this threshold, 
provides a simple and accurate way of identifying doublets. 

A second aspect of the present invention provides a method of seismic data acquisition 
comprising: acquiring first seismic data at a first seismic receiver and simultaneously 
acquiring second seismic data at a second seismic receiver spatially separated from the 
first seismic receiver, and processing the first and second seismic data according to a 
method of the first aspect of the invention. 

A third aspect of the invention provides an apparatus for processing seismic data, the 
seismic data comprising first data acquired at a first seismic receiver and second seismic 
data acquired at a second receiver spatially separated from the first seismic receiver, the 
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apparatus comprising means for determining an overall measure of similarity for a pair 
of events in the seismic traces, the overall measure of similarity being indicative of 
similarity between the events within the first data traces and of similarity between the 
events within the second data traces. 

5 

The apparatus may comprise a programmable data processor. 

A fourth aspect of the present invention provides a storage medium containing a 
program for the data processor of an apparatus according to the third aspect of the 
10 invention. 

A fifth aspect of the invention provides a storage medium containing a program for 
controlling a programmable data processor to carry out a method of the first aspect of 
the invention. 
15 " 

A sixth aspect of the invention provides a program for controlling a computer 
programmable data processor to carry out a method of the first aspect of the invention. 

Further preferred features of the invention are set out in the dependent claims. 

20 

Brief Description of the Drawings 

Preferred embodiments of the present invention will now be described by way of 
illustrative example with reference to the accompanying figures, in which: 

25 

Figure 1 is a schematic illustration of the occurrence of doublets in seismic data; 
Figure 2 is a schematic flow diagram of a method of the invention; 
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Figure 3(a) is a flow diagram of one method of determining the overall cross-correlation 
coefficient; 



r 8 

Figure 3(b) illustrates one step of the method of Figure 3(a); 

Figure 4 illustrates results obtained by a method of the present invention when applied 
to a synthetic dataset; and 

5 

Figure 5 is a block schematic diagram of a data processing apparatus according to the 
present invention. 

Detailed Description of the Preferred Embodiments- 

10 

A preferred embodiment of the invention will be described in which the overall measure 
of similarity is an overall correlation coefficient. However, the invention is not limited 
to use of an overall correlation coefficient as the overall measure of similarity and, in 
principle any quantity which gives a measure of the similarity between two waveforms 
15 may be used. 

This embodiment of the invention is based on waveform correlation. The cross- 
correlation coefficient c xy between a first waveform x(t) and a second waveform y(t) is a 
• dimensionless quantity that gives a measure of the similarity between the waveform x(i) 
20 and the waveform y(t). The cross-correlation coefficient is defined as: 

oo 

$x(t)y(t + r)dt 

c^T) = j= (1) 

In equation (1), x(t) is the time series denoting the first waveform, and y(t) is the time 
25 series denoting the second waveform. y(t + r) denotes the time series for the second 
waveform shifted by a delay time r. 

The cross-correlation coefficient is a function of the delay time r. If x(t) and y(t) are 
generally similar waveforms, but one is time-shifted relative to the other by a time To, 
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the cross-coirelation coefficient of the two waveforms would contain a maximum at a 
delay time of t = x 0 . The magnitude of the cross-correlation coefficient takes a value 
between 1 and -1, with a" value of 1 indicating exact correlation between the two 
waveforms and a value of zero indicating that the two waveforms are completely un- 
5 correlated with one another. A negative value of the cross-correlation, coefficient 
indicates that the two waveforms are anti-correlated (i.e., that one waveform is 
correlated with the negative of the other waveform). 

In equation (1) above the cross-correlation coefficient is calculated in the time domain. 
10 It may alternatively be calculated in the frequency domain, by analogy with the 
convolution theorem, according to: 



c xy (.T)=±F D - l [X\f)nf)] (2) 

15 In equation (2), F D ' 1 denotes the inverse discrete Fourier transform, X*(f) is the complex 
conjugate of the Fourier transform of x(t), Y(f) is the Fourier transform of y(t) and N is 
the number of data points. 

* . 
i. 

Figure 2 is a schematic diagram of one method of the present invention. 

20 

Initially at step 1, seismic data traces are acquired at at least two spatially separated 
seismic receivers. An arrangement of three receivers, as shown in Figure 1, may be 
used but the invention is not limited to three receivers and may alternatively be carried 
out using two receivers or using more than three receivers. At least one seismic data 

25 trace is acquired for each receiver. In general, each receiver will acquire more than one 
trace, with each trace representing a different component of the seismic wavefield - and 
commonly each receiver will acquire three traces, representing the x-, y- and z- 
components. The seismic data traces are acquired simultaneously at each receiver, so 
that the emission of acoustic energy in the survey area will, in general, produce an event 

30 in the seismic data traces acquired at each of the receivers. 
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The method may alternatively be applied to pre-existing seismic data. In this case, step 
1 is not carried out and is replaced by step 2 of retrieving pre-existing seismic data from 
storage. 

5 At step 3, the overall cross-correlation coefficient of the acquired seismic data, for a pair 
of seismic events, is determined. This step is described in more detail below. As 
indicated above, the overall cross-correlation coefficient is indicative of the correlation 
between a pair of events within the data traces acquired at each receiver. 

10 At step 4, the overall cross-correlation coefficient for the pair of seismic events is 
compared with a pre-determined threshold. The threshold may be defined by the user 
and may depend on, for example, the noise in the seismic data traces. For example, a 
user may decide that, for a particular set of acquired seismic data, an overall cross- 
correlation coefficient of 0.9 or greater indicates a doublet and, in this case, step 4 

15 comprises determining whether the overall cross-correlation coefficient is equal to or 
greater than 0.9. If the result of the step 4 is that the overall cross-correlation coefficient 
is less than the pre-determined threshold, it is judged that the pair of events do not form 
a doublet. A "y es " determination at step 4 may be taken as indicating that the pair of 
seismic events form a doublet. 

20 

In this embodiment, therefore, a pair of events is identified as a doublet if the overall 
correlation coefficient is equal to or greater than a pre-set threshold. As noted above, 
the method of the invention may be performed in real-time or near real-time, and the 
thus invention makes possible real-time identification of doublets. 

25 

The method of Figure 2 may be repeated for each pair of events in the seismic data 
traces. 

One method, of determining the overall cross-correlation coefficient of the invention 
30 will now be described, with reference to Figure 3(a). The method will be described - 
with particular reference to the case where the seismic data are or were acquired using 
two "or more spatially separated multi-component seismic receivers. A multi-component 
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seismic receiver contains, as is well-known, two or more sensors, with each sensor 
arranged to record a different component of the seismic wavefield incident on the 
receiver. In a land-based survey a 3-component seismic receiver is commonly used, and 
this is provided with three mutually orthogonal sensors. In general, a 3-component 
receiver is oriented in use so as to record the waveform components in the x-, y- and z- 
directions (with the z-direction being defined as the vertical direction). Thus, it is 
assumed that each of the receivers 7, 8, 9 is able to record three data traces, DJt), D yi (t) 
and Ddt) which are respectively the x-, y- and z-components of a parameter of the 
seismic wavefield (such as the x-, y- and z- components of the particle velocity) 
acquired at the I th receiver. Each of these data traces extends from t = 0 to t = T, and 
records the amplitude and phase of the respective component as a function of time from 
t = 0 to t = T. 

The first step in processing the seismic data is to compute cross-correlation coefficients 
for each of the seismic components acquired at the I th receiver. Thus, at step 1 of Figure 
3(a) the trace D„{t), the x-component data for the i* receiver, is investigated. The trace 
will contain a number of events as shown schematically in Figure 3(b), and at step 1 the 
cross-correlation coefficient between a time window 2} containing event; and a time 
window T k containing event k is determined. This produces the cross-correlation 
coefficient C^Jt) - me cross-correlation coefficient between events j and k in the x- 
component seismic data trace acquired at the i* receiver. The cross-correlation 
coefficient d k Jr) is a function of time lag and extends over the time-lag range from x = 
0 to x = T. In general, C^r) will have a maximum or minimum for a certain value of 
the time lag. The time lag at which the cross-correlation coefficient has a 

maximum or minimum indicates the time delay between event j and event k, and the 
value of the cross-correlation coefficient at its maximum or rmnimum is indicative of 
the degree of correlation between event; and event k. If events jjc form a doublet they 
would have similar waveforms and so the maximum value of the cross-correlation 
coefficient C^C ^ would be near to 1. 

The cross-correlation coefficient d^r) for the events ; and k in D^r) - the y- 
component seismic data trace for the I th receiver - is similarly obtained by cross- 
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correlating the time window containing event j in the y-component data trace with the 
time window containing event k in the y-component data trace at step 2. The cross- 
correlation coefficient d k ^T) for the events j and k in the z-component seismic data 
trace for the I th receiver - DJrf - is obtained by cross-correlating the time window 
5 containing event j in the z-component data trace with the time window containing event 
k in the zrcomponent data trace at step 3. Steps 1, 2 and 3 do not need to be carried out 
in the order described, and could in principle be performed simultaneously. 

The cross-correlation coefficients d k ^r), C jk yi (T), d^fr) are preferably calculated in 
10 the frequency domain, according to equation (2) since, for a given computer, this is 
generally faster than calculating the cross-correlation coefficients in the time domain. 
In principle, however, the cross-correlation coefficients could be calculated in the time 
domain. 

15 Next, at step 4, a combined cross-correlation coefficient for the events j and k in the data 
traces * acquired by the I th receiver is derived from the individual cross-correlation 
coefficients for these events in the x-, y- and z-components of the seismic data acquired 
by the I th receiver. In one embodiment, the combined cross-correlation coefficient for 
the events j and k for the I th receiver is simply obtained by averaging the individual 

20 cross-correlation coefficients obtained for these events in the x-, y- and z- components 
of the seismic data acquired at the i a receiver. In a preferred embodiment, however, a 
weighted averaging process is used, in which the combined cross-correlation coefficient 
for the events j and k for the I th receiver, d k Rl {t) y is determined according to: 



In equation (3), d k ^T), d k yi {T) and d k zl {r) are,- as described above, the cross-correlation 
coefficients for events j and k in, respectively, the x- y y- and z-components of the 
seismic data acquired at the i* receiver, and A x , A y and A z are the maximum amplitudes 
30 for, respectively, the jc-, y- and z-components of the seismic data acquired at the z ,th 
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C£(r) = 



A X C* (r) + A y Cff (r) + A Z C* (r) 
A x + A y + A z 



(3) 
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receiver. The cross-correlation coefficient C*« is again a function of time-lag and 
extends over the time-lag range from x = 0 to x = T. 

The effect of the averaging step is to increase the weight given to a pair of events that 
occur, with the same time-lag or a similar time-lag, in the x-, y- and z-components of the 
data acquired at the i* receiver. A pair of events having a large cross-correlation at a 
certain time-lag in one component of the data, but which do not have a large cross- 
correlation coefficient in other components of the data, will have its importance reduced 
by the averaging process. Use of the weighted average of equation (3) reduces the 
effect of noise on the determination of the cross-correlation coefficient d k Rl {r), because 
of the amplitude weighting introduced by equation (3). 

Next, it is determined whether the cross-correlation coefficient for the events; and k has 
been calculated for all receivers, at step 5. If step 5 yields a "no" determination, a 
suitable counter is incremented at step 6, and steps 1, 2, 3 and 4 are repeated to 
determine the cross-correlation coefficient C jk R , M (T)fox the events / and k for the (i+lf 
receiver. Steps 5, 6 and 1 to 4 are repeated until a cross-correlation coefficient for the 
events j and k has been determined for each receiver. 

The next step is to determine an overall cross-correlation coefficient from the cross- 
correlation coefficients for each receiver. As noted above in a seismic surveying 
arrangement having closely spaced receivers (for example, a receiver separation of 
around 50m) the two events of a doublet should occur in seismic data acquired at all 
receivers with substantially the same time-lag for all receivers - and a pair of events that 
are recorded with a similar time-lag at each receiver will contribute to a high value for 
the overall cross-correlation coefficient. However, a pair of events occurring with a 
certain time-lag in the data acquired at one receiver but for which there is no 
corresponding pair of events at a similar time-lag in the data acquired at other receivers 
will not produce a high value for the overall cross-correlation coefficient. 

In one embodiment of the invention, the overall cross-correlation coefficient d k (r) for 
the events j and k is obtained by finding the average value of the coefficients for 
each receiver. The time-lag between the two events of a doublet is unlikely to be 
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exactly the same in the data acquired by all receivers, however, since the source- 
receiver distance will be different for different receivers. Simply summing the value of 
the individual d k Ms will generally not be effective therefore, since the delay time at 
which d k Ri (r) is a maximum will vary slightly from one receiver to another. In a 
preferred embodiment, therefore, the averaging process is carried out within a time 
window which takes account of the expected variation in time lag from one receiver to 
another. This is shown schematically in Figure 3(c). The time window may be defined 
by a user, taking account of, for example, the receiver separation and the velocity of 
seismic energy in the reservoir layer. 

In this embodiment, the overall cross-correlation coefficient C? k (x) for the events j and k 
is obtained by finding the delay time at which the average, in the chosen time window, 
of the coefficients d k Ri (r) for each receiver is greatest. That is: 
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C*(r) = - ~ 



max 



/=1 



(4) 



m 



In equation (4) d k Ri is the cross-correlation coefficient for the events j and k for the I th 
receiver, m denotes the total number of receivers, and %i e fa-St, xH-5t), for I = 1,. . .m. 

20 Equation (4) includes a time window of 28t for the time lag which, as noted above, may 
be defined by the user. A pair of events that occur within this window for the time-lag, 
in the data acquired at all receivers, will contribute to a maximum in the overall cross- 
correlation coefficient C(t). 

25 Determination of the overall cross-correlation coefficient C Jk (r) is step 7 in Figure 3(a). 



This process is then repeated for other pairs of different events in the seismic data 
traces. (The correlation coefficient of a waveform with itself always has a peak value of 
1 at zero time lag. The method is therefore preferably arranged to avoid determining the 
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correlation coefficient of an event with itself at steps 1, 2 and 3.) This is represented 
schematically in Figure 3(a) by step 8 of determining whether an overall cross- 
correlation coefficient has been determined for each pair of different events in the data 
traces and step 9 of selecting another pair of different events (such as event k and event 
Z) in the case of a "no" determination at step 8. Steps 1-7 are then repeated for the new 
pair of events. 

Figure 4 illustrates results of the method of the present invention. The results of figure 
4 were obtained by applying the above method to a synthetic data set synthesised for a 
range of known event locations and a range of known source mechanisms. The data 
was synthesised for three receivers, arranged along a vertical line to simulate a 
borehole. 

The simulation used 20 sources arranged at different spatial locations. The x-, y- and z- 
components of the seismic data acquired at each receiver as a result of emission of 
seismic energy at each source location were simulated. 

The simulated seismic data were then cross-correlated. The cross-correlation was 
carried out according to the method described above - that is, a respective cross- 
correlation coefficient was determined for a pair of events in each of the x-, y- and Z- 
components of the seismic data simulated for one receiver, and the combined cross- 
correlation coefficient for the events, for that receiver, was obtained as a weighted 
average of the respective cross-correlation coefficients using equation (3). The 
combined cross-correlation coefficient for the second receiver and third receivers were 
obtained in a similar way, and the overall cross-correlation coefficient for the three 
receivers was obtained using equation (4). 

This process was repeated for every combination of two sets of data simulated for 
different source locations. 



Figure 4 shows the greatest magnitude of the overall cross-correlation coefficient, 
C jk (x), plotted against the event separation. (The event separation is the separation 
between the locations of the two sources used in the cross-correlation). 

5 The two events of a doublet would be expected to have an event separation of no more 
than one quarter of the dominant seismic wavelength. In the simulated data used to 
obtain figure 4, the dominant frequency was approximately 100Hz and, assuming a 
maximum seismic velocity of 6000m/s, one quarter of the dominant wavelength is 
approximately 15m. It is therefore expected that doublets should have an event 
10 separation of less than approximately 15m. 

Figure 4 shows that, at event -separations of less than 15m where all of the events are 
assumed to be doublets, the overall cross-correlation coefficient C(t)s generally has a 
maximum value of greater than 0.9. At event separations of greater than 15m, where 
15 the events are assumed to be non-doublets, the overall cross-correlation coefficients are 
generally lower than 0.9. It was also found that, at event separations greater than 15m, 
the average of the overall cross-correlation coefficients is approximately zero indicating 
that, on average, no correlation exists. 

20 It may be seen from figure 4 that one criteria suitable for distinguishing doublets from 
non-doublets would be to determine whether the overall cross-correlation coefficient is 
greater than 0.9. If the overall cross-correlation coefficient is > 0.9, it can be assumed 
that the pair of events fonn a doublet. In the results of figure 4, this criterion would 
correctly identify 97% of the doublets. 

25 

It will be seen in figure 4 that the criterion of an overall cross-correlation coefficient of 
greater than 0.9 identifies a few event pairs that have event separations of between 15 
and 50m as doublets. Since these events have an event separation greater than one 
quarter of the dominant wavelength they would not generally be regarded as doublets. 
30 They have travelled through similar earth structure, however, since the synthetic model 
is much simpler than the real structure of the earth. Therefore, the events may be 
considered as doublets for the simple earth model used in the calculations. 
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The results of Figure 4 were obtained using noise-free synthetic seismic data traces. 
The effects of noise- were investigated by adding noise to the synthetic data traces. It 
was found that, with increasing noise level, there was a general reduction in the 
magnitude of the overall cross-correlation coefficient C(x) and that it became harder to 
distinguish doublets from non-doublets. For 5% added noise, results similar to those of 
Figure 4 were obtained, and the 0.9 threshold for identifying muliplets was -still 
applicable. For 10% added noise, the threshold for identifying doublets had dropped to 
0.85, and increasing numbers of events With separations than 15m were identified as 
doublets. When a noise level of around 20% was used, the threshold for identifying 
doublets had dropped further to 0.65. With a 50% noise level it was not possible to 
define a threshold for identifying doublets, even though there was still a noticeable 
difference in the nature of the cross-correlation coefficient for doublets compared to 
non-doublets. It was difficult to apply the method of the invention for noise levels 
significantiy greater than 20%. 1 ' 

Figure 5 is a schematic block diagram of a programmable apparatus 10 according to the 
present invention. The apparatus comprises a programmable data process 11 with a 
programme memory 12, for instance in the form of a read-only memory (ROM), storing 
a programme for controlling the data processor 11 to perform any of the processing 
methods described above. The apparatus further comprises non-volatile read/write 
memory 13 for storing, for example, any data which must be retained in the absence of 
power supply. A "working" or scratch pad memory for the data processor is provided 
by a random access memory (RAM) 14. An input interface 15 is provided, for instance 
for receiving commands and data. An output interface 16 is provided, for instance for 
displaying information relating to the progress and result of the method. Seismic data 
for processing may be supplied via the input interface 15, or may alternatively be 
retrieved from a machine-readable data store 17. 



The programme for operating the system and for performing the method described 
hereinbefore is stored in the programme memory 12, which may be embodied as a semi- 
conductor memory, for instance of the well-known ROM type. However, the 
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programme may be stored in any other suitable storage medium, such as magnetic data 
carrier 12a, such as a "floppy disk" or CD-ROM 12b. 
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CLAIMS: 

1. A method of processing seismic data, the seismic data comprising at least first 
seismic data traces acquired at a first seismic receiver and second seismic data traces 
acquired at a second receiver spatially separated from the first receiver, the method 
comprising: determining an overall measure of similarity for a pair of events in the 
seismic traces, the overall measure of similarity being indicative of similarity between 
the events acquired at the first seismic receiver and of similarity between the events 
acquired at the second seismic receiver. 

2. A method as claimed in claim 1 wherein the method is a real-time processing 
method. 

3. A method as claimed in claim 1 wherein the overall measure 6f similarity is an 
overall correlation coefficient indicative of correlation between the events acquired at 
the first seismic receiver and of correlation between the events acquired at the second 
seismic receiver. 

4. A method as claimed in claim 3 wherein determining the overall cross- 
correlation coefficient comprises 

(a) determining a first correlation coefficient for the pair of events from the first 
data; 

(b) determining a second correlation coefficient for the pair of events from the 
second data; and 

(c) determining the overall correlation coefficient for the pair of events from the 
first correlation coefficient and the second correlation coefficient. 

5. A method as claimed in claim 1, 2, 3 or 4 and comprising the further step of 
comparing the overall measure of similarity for the pair of events with a first pre- 
determined threshold. 
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6. A method as claimed in claim 5 and comprising the further step of identifying 
the pair of seismic events as a doublet if the overall measure of similarity is equal to or 
greater then the. first threshold 

5 7. A method as claimed in claim 4, or in claim 5 or 6 when directly or indirectly 
dependent from claim 4, wherein the first seismic receiver is a multi-component seismic 
receiver, and step (a) comprises: 

(i) determining respective correlation coefficients for the pair of events for each 
data component acquired by the first seismic receiver; and 
10 (ii) determining the first correlation coefficient for the pair of events from the 
respective correlation coefficients. 

8. A method as claimed in claim 7 wherein step (i) comprises determining the 
respective correlation coefficients in the frequency domain. 

15 

9. A method as claimed in claim 7 or 8 wherein step (ii) comprises determining the 
first correlation coefficient as a weighted average of the respective correlation 
coefficients. 

20 10. A method as claimed in claim 4, or in claim 5 or 6 when directly or indirectly 
dependent from claim 4, wherein the second seismic receiver is a multi-component 
seismic receiver, and step (b) comprises: 

(i) determining respective correlation coefficients for the pair of events for each 
data component acquired by the second seismic receiver; and 
25 (ii) determining the second correlation coefficient for the pair of events from the 
respective correlation coefficients for each data component acquired by the second 
seismic receiver. 

11. A method as claimed in claim 10 wherein step (i) comprises determining the 
30 respective correlation coefficients in the frequency domain. 
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12. A method as claimed in claim 10 or 11 wherein step (ii) comprises determining 
the second correlation coefficient as a weighted average of the respective correlation 
coefficients for each data component acquired by the second seismic receiver. 

13. A method as claimed in claim 4, or in claim 5 or 6 when directly or indirectly 
dependent from claim 4, wherein step (c) comprises determining the overall correlation 
coefficient for the pair of events as an average of the first correlation coefficient for the 
pair of events and the second correlation coefficient for the pair of events. 

14. A method as claimed in claim 13 and comprising determining the overall 
correlation coefficient for the pair of events according to: 



max 

C(r) = 



m 



where Cw is the cross-correlation coefficient for the pair of events for the i* receiver, 
t, e (r, —ki,t l + At) , for I = 1,. . .m and m is the number of receivers. 

15. A method of seismic data acquisition comprising: acquiring first seismic data at 
a first seismic receiver and simultaneously acquiring second seismic data at a second 
seismic receiver spatially separated from the first seismic receiver; and processing the 
first and second seismic data according to a method as defined in any of claim 1 to 14. 

16. An apparatus for processing seismic data, the seismic data comprising first data 
acquired at a first seismic receiver and second seismic data acquired at a second seismic 
receiver spatially separated from the first receiver, the apparatus comprising means for 
determining an overall measure of similarity for a pair of events in the seismic traces, 
the. overall measure of similarity being indicative of similarity between the events 
within the first data traces and of similarity between the events within the second data 
traces. 
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17. An apparatus as claimed in claim 16 and adapted to process the seismic data in 
real time. 

18. An apparatus as claimed in claim 16 or 17 and adapted to determine, as the 
5 overall measure of similarity, an overall correlation coefficient indicative of correlation 

between the events within the first data traces and of correlation between the events 
within the second data traces. 

19. An apparatus as claimed in claim 18 and comprising: 

10 means for determining a first correlation coefficient for the pair of events from 

the first data; 

means for determining a second correlation coefficient for the pair of events 
from the second data; 

means for determining the overall correlation coefficient for the pair of events 
15 from the first correlation coefficient and the second correlation coefficient. 

20. An apparatus as claimed in any of claims 16 to 19 and further comprising means 
for comparing the overall measure of similarity for the pair of events with a first pre- 
determined threshold. 

20 

21. An apparatus as claimed in claim 16, 17, 18, 19 or 20 and comprising a 
programmable data processor. 

22. A storage medium containing a program for the data processor of an apparatus 
25 as defined in claim 21 . 

23. A storage medium containing a program for controlling a programmable data 
processor to carry out a method as defined in any of claims 1 to 14. 

30 24. A program for controlling a computer to carry- out a method as defined in any of 
claims 1 to 14. 
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ABSTRACT 
Processing Seismic Data 



A method of processing seismic data that contains at least first seismic data traces 
acquired at a first seismic receiver and second seismic data traces acquired at a second 
receiver spatially separated from the first receiver comprises determining an overall 
measure of similarity for a pair of events in the seismic traces. The overall measure of 
similarity is indicative of similarity between the events acquired at the first seismic 
receiver and of similarity between the events acquired at the second seismic receiver. 



In one method, the overall measure of similarity is an overall cross-correlation 
coefficient. The oVerall cross-correlation coefficient is found by determining a first 
correlation coefficient for the pair of events from the data acquired at the first receiver 
and determining a second correlation coefficient for the pair of events from the data 
15 acquired at the second receiver. The overall correlation coefficient for the pair of events 
may be obtained from the first correlation coefficient and the second correlation 
coefficient by an averaging process. 

The overall measure of similarity may be compared with a threshold to determine 
20 whether the pair of events form a doublet. 

The method makes possible real-time or near-real-time identification of doublets. 



[Figs 2, 3(a)] 

25 




7— to 
■z. — a 



w 

c 



M&C Folio No P52475GB 



acquire data at at 
least two receivers 



determine overall 
correlation coefficient C(x) 



retrieve data 
from storage 




not Multiplet 



Multiplet 



M&C Folio No P52475GB 



23 



Cx(t) for I th receiver 



Cy(t) for ^receiver 



Cz(t) for ^receiver 



C> i-n+1 



C R for i* 1 receiver j 



select new 
event pair 






i 


determine C(t) 




r 




M&C Folio P52475GB . , *> - 

c 




threshold = 0.9 



+ 
+++ 

X ++ 

> ++ 

+ + + 

+ + 
+ 

+ 



+ 



J + * ++ 
* ++ 



+ 

+ + 



3D 40 50 60 
Event Separation 



BO 



9D 1D0 



Figure 4- 



10 



]~Data ~! 
' store j" 



is 



il 



Input 




interface 





2l 



Programmable 
data processor 



"Workmg" 




memory 




[RAM] 





Non -volatile 
read/write 
memory 



1 



\1o 



Output interface 





Program 




memory 




[ROM] 



«2 5 ; 



'Floppy] 
1 disk 1 



,J ^/-iCD-ROM| 



FIG 5" 



PCT/GB2004/002614 



This Page is Inserted by IFW Indexing and Scanning 
Operations and is not part of the Official Record 

BEST AVAILABLE IMAGES 

Defective images within this document are accurate representations of the original 
documents submitted by the applicant. 

Defects in the images include but are not limited to the items checked: 

□^BfcACK BORDERS 

LJ IMAGE CUT OFF AT TOP, BOTTOM OR SIDES 

□ FADED TEXT OR DRAWING 
□ / BLURRED OR ILLEGIBLE TEXT OR DRAWING 
jZfsKEWED/SLANTED IMAGES 

□ COLOR OR BLACK AND WHITE PHOTOGRAPHS 

□ GRAY SCALE DOCUMENTS 

□ LINES OR MARKS ON ORIGINAL DOCUMENT 

□ REFERENCE(S) OR EXHIBIT(S) SUBMITTED ARE POOR QUALITY 

□ OTHER: 

IMAGES ARE BEST AVAILABLE COPY. 
As rescanning these documents will not correct the image 
problems checked, please do not report these problems to 
the IFW Image Problem Mailbox. 



